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Extensive density-functional calculations are performed to understand atomic chemisorption on 
the TiC(lll) and TiN(lll) surfaces, in particular the calculated pyramid-shaped trends in the ad- 
sorption energies for second- and third-period adatoms. Our previously proposed concerted-coupling 
model for chemisorption on TiC(lll) is tested against new results for adsorption on TiN(lll) and 
found to apply on this surface as well, thus reflecting both similarities and differences in electronic 
structure between the two compounds. 
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I. INTRODUCTION 

Within an extensive theoretical study of atomic 
chemisorption on the Ti-terminated polar TiX(lll) (X 
= C, N) surfaces, we have presented characteristic trends 
in adsorption energies and geometries for H, second-row 
elements B, C, N, O, and F, and third-row elements Al, 
Si, P, S, CI P, 0]- The study has both fundamental and 
technological motivations, as these materials are impor- 
tant in applications in which surface properties are essen- 
tial, such as growth of alumina on Ti(C,N)(lll) and use 
of TiN as a thin diffusion barrier in integrated circuits or 
as a biocompatible material. 

Our results show similar trends for TiC and TiN, with 
atomic adsorption energies for the second- and third- 
period elements varying along the periods in a charac- 
teristic pyramid- like way [Fig. |3Ia)] Within each 
period, the energies increase from group III to group VI 
and decrease from group VI to group VII, with over- 
all higher values for period 2. These variations are very 
large, from ~ 3.5 eV for Al to ~ 8.5 eV for O. Also, the 
preference for adsorption on TiC vs. TiN varies, with a 
slightly stronger chemisorption on TiN for the elements 
toward the left of the periodic table (B and Al) and a 
slightly stronger chemisorption on TiC for the remaining 
adatoms. 

Such variations cannot be explained within a simple 
model, where only one type of coupling mechanism is 
considered. In a separate study, we propose a concerted- 
coupling model for adsorption on TiC(lll), in which two 
different types of adatom-substrate couplings are active 
. This is similar to the c?-band model for chemisorp- 
tion on transition metals 3], however, with very different 
couplings, due to the more complex electronic structure 
of TiC. In this contribution, we pursue our investigation 
further, analyzing the electronic structure of the different 
adsorbates on TiN(lll) and discussing whether our pro- 
posed model can be generalized and used to understand 
the similarities and differences between the TiN(lll) and 
TiC(lll) substrates. 
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II. COMPUTATIONAL METHOD 

The computational method used is described in Ref. 
2]. In short, the pseudopotential-plane-wave density- 
functional-theory (DFT) code dacapo is used, with 
PW91 GGA exchange-correlation potential, ultrasoft 
pseudopotentials, and slab geometry 

In the present paper, the electronic structures of the 
considered systems are analyzed by calculating and plot- 
ting the total and local densities of states (DOS and 
LDOS, respectively) obtained from the calculated Kohn- 
Sham (KS) wavefunctions and energy eigenvalues. Atom- 
and/or orbital-projected LDOS(i?)'s are obtained by pro- 
jecting the KS wavefunctions onto the individual atomic 
orbitals and then plotting them as functions of the energy 
relative to the Fermi level, E — Ep. The spatial extent of 
the DOS's at specific energy intervals is also analyzed in 
detail by examining the KS-wavefunction amplitudes for 
the relevant electronic states in three-dimensional real 
space. Also, the charge localization around individual 
atoms is analyzed with the "atoms-in-molecule" method 
of Bader Q. 

III. RESULTS AND ANALYSIS 

A. Electronic structure of the clean TiX(lll) 
surfaces 

Bulk TiC and TiN adopt the NaCl structure. Thus, in 
the (111) direction, the crystals are composed of alter- 
nating Ti and X atomic layers, in an ABC stacking se- 
quence. Our study focuses on the Ti-terminated surfaces, 
since experiments show this to be the stable termination 
for TiC(lll) 0. The calculated lattice parameters and 
surface relaxations are given in Refs. P, 

Both Ti-terminated TiX(lll) surfaces are character- 
ized by the presence of strong Ti-localized surface reso- 
nances (TiSR's) (7(1, pinned at Ep [Figs.QIa-b)] and con- 
sisting of mainly d(^xz,yz) ^-^d '^(ij/.s^-yS) orbitals (where 
z is perpendicular to the surface). Figure ^c) shows this 
symmetry and how the TiSR electron distribution is lo- 
calized mainly around the Ti atoms (as evidenced also by 
the small width of the TiSR peaks) and extends toward 
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the fee surfaee sites. 

Our ealeulated surface DOS's [Figs.^a-b)] show that 
the TiSR is more filled on TiN(lll) than on TiC(lll). 
This is eonfirmed by a Bader analysis of the surface 
charges, yielding a surface-Ti ionicity of +1.09 for TiC 
and of +0.96 for TiN. This can be understood by con- 
sidering the filling of the TiSR's as the mechanism that 
compensates the electrostatic instability of the polar 
TiX(lll) surfaces 0. As bulk TiN is more ionic than 
bulk TiC (Ti — > X charge transfers of 1.51 electrons for 
TiC and 1.62 for TiN, according to our Bader analysis) 
0, a larger compensating surface charge is needed on 
TiN(lll). 

The bulk electronic structure of TiX [inserts in Figs, 
n^a-b)] consists of (i) a lower valence band (LVB), of only 
X2s character, (ii) an upper valence band (UVB), char- 
acterized by bonding Ti3c?-X2p states in its high-energy 
part and by X2p-X2p bonding interactions in its low- 
energy part, and (iii) a conduction band (CB), of mainly 
Ti3d character, corresponding to antibonding Ti3d-X2p 
states P, • The ionicity of bulk TiX is evidenced by a 
predominance of Ti (X) character in the CB (UVB). 

Both TiX(lll)-surface UVB's [Figs. IHa-b)] are com- 
posed of one main peak, three smaller peaks at lower 
energies, a shoulder on the high-energy side of the main 
peak, and a small peak at the upper edge of the UVB. 
For both surfaces, our detailed analysis of the individual 
KS wavefunctions composing the UVB (Fig. ^ shows 
that, similarly to the bulk UVB, the main UVB peak 
and the high-energy UVB region consist of Ti-X bond- 
ing states [Fig. [2Ic)]. On the other hand, the three low- 
energy peaks consist almost exclusively of X-X bond- 
ing states that are mainly localized to the surface region 
[Figs. Efa-b)], that is, they are X-localized surface reso- 
nances (XSR's). 

Like for the bulk DOS 0, several differences exist be- 
tween the TiN(lll) and TiC(lll) surface DOS's [Figs. 
IHa-b)]: (i) the UVB and CB of TiN(lll) He at lower 
energies, due to the partial filling of the CB with the 
extra electron per formula unit present in TiN; (ii) the 
pseudogap between the UVB and the TiSR is larger for 
TiN(lll); and (iii) the UVB DOS is considerably larger 
and more X-localized on TiN, which again is an indica- 
tion of the higher ionicity of TiN. 

Thus, both surface DOS's differ from their correspond- 
ing bulk DOS's [see Figs. CJa-b)] by (i) the presence of 
TiSR's around Ep, (ii) a quenching of the UVB peak ly- 
ing above the main peak, and (iii) the presence of XSR's 
in the low-energy UVB region. These changes can be 
understood in terms of breakage of the Ti-X and X-X 
bonds that cross the (111) plane upon cleavage of the 
bulk structures to create the (111) surfaces. As a conse- 
quence, the bonding and antibonding Ti-X states, which 
are located in the high-energy part of the UVB and in 
the CB, respectively, collapse into more atomic-like X 
and Ti orbitals. Therefore, the quenching of the high- 
energy UVB peak can be associated with the disappear- 
ance of X-localized bonding Ti-X states, while the ap- 



pearance of the TiSR's arises from the energetical lower- 
ing of Ti-localized antibonding Ti-X states (indeed, for 
both compounds, a closer inspection of the CB reveals 
the quenching of a bulk CB peak located right above the 
TiSR's). The TiSR's can thus be interpreted as dangling 
bonds extending toward the fee sites (which in the bulk 
would be occupied by X atoms) [Fig.^c)]. Similarly, the 
appearance of the XSR's in the low-energy UVB region 
is caused by the breakage of X-X bonding interactions, 
which dominate at these energies. 



B. Electronic structure of fee adatoms on TiX(lll) 

In a separate study, we have proposed a model for the 
description of the exceptionally strong atomic chemisorp- 
tion on TiC(lll) Q. In this concerted-coupling model, 
based on the Newns- Anderson (NA) model the ad- 
sorption is understood to arise from the concerted action 
of two different types of adatom-substrate coupling: a 
strong adatom-TiSR coupling and weaker adatom-CSR 
couplings. 

Briefly, the NA model for chemisorption on metal 
surfaces yields two limiting cases: weak and strong 
chemisorption, corresponding, respectively, to small and 
large coupling matrix elements between adatom and sub- 
strate states, compared to the substrate band width. In 
the weak chemisorption, the adsorbate level is shifted 
and broadened (if it lies in the energetical region of 
the substrate band). In the strong chemisorption, the 
adatom level is replaced by two levels that lie on each 
side of the substrate band. These are usually interpreted 
as the bonding and antibonding states of the "surface 
molecule" formed by the adatom and its neighboring sub- 
strate atoms. These bonding and antibonding states ap- 
pear gradually, as the adatom-substrate hopping integral 
h increases. At small h, the adlevel is broadened. As h 
increases, this level broadens further and splits gradually 
into two separate levels TO]. 

Our concerted-coupling model for TiC(lll) is based 
on an analysis of the difference in DOS (ADOS) be- 
fore and after adsorption (thus, the ADOS shows the 
appearance and quenching of states due to the adsorp- 
tion of an adatom as positive and negative peaks, re- 
spectively) 0. For all considered adatoms on TiC(lll), 
our ADOS's show: (i) the appearance of Ti-dominated 
states at ~ +1.0 eV; (ii) a strong quenching of the TiSR; 
(iii) the appearance of states below Ep, whose energies 
decrease with increasing adatom number Z within each 
adatom period (due to the decreasing free-atom energy 
as Z increases). 

Within the NA model, these results are interpreted as 
evidence for a strong adatom-TiSR coupling, giving rise 
to antibonding and bonding states above and below, re- 
spectively, Ep. At the same time, our calculated ADOS's 
show also that the bonding adatom-TiSR level below Ep 
is successively broadened and split into subpeaks as its 
energy approaches the middle of the substrate UVB re- 
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FIG. 1: (a-b) Total and atom-projected DOS for the top TiX bilayer of the clean (a) TiC(lll) and (b) TiN(lll) surfaces. The 
inserts in (a) and (b) show the total DOS for bulk TiC and TIN, respectively, (c) Real-space contour plot (from above the 
surface and showing only the top TiC bilayer) of the total charge density of the electronic states above —0.6 eV on TiC(lll), 
showing the TiSR (C atoms: green balls; Ti atoms: inside the electron clouds, as marked in the figure) . 




FIG. 2: Real-space contour plots of the KS wavefunctions for the clean TiN(lll) surface at selected energies, viewed perpen- 
dicularly to the surface, illustrating: (a-b) typical bonding N-N interactions, with and without coupling to bulk states, in the 
lower UVB region; (c) a typical Ti-N bonding state, with high N-localization, in the main UVB peak. Ti (N) atoms are grey 
(green) balls, or lie inside the electron clouds, as marked in the figures. 



gion. This indicates, within the NA model, an additional, 
weaker, coupling between the TiSR-modified adlevel and 
the substrate UVB peaks, which are characterized by 
CSR's. Thus, the subpeaks correspond to combinations 
of the bonding and antibonding solutions that arise from 
adatom-CSR couplings. This is confirmed by (i) the 
presence, in the substrate-projected ADOS's, of nega- 
tive peaks that are almost exclusively localized around 
the C atoms and that lie at energies in between the sub- 
peaks and (ii) our state-resolved, real-space, analysis of 
the DOS, which shows that the subpeaks are dominated 
by adatom-C bonding states at lower UVB energies and 
by adatom-Ti bonding states at higher energies. The 
adatom-C states are only present at energies lower than 
the main substrate UVB peak and can only be inter- 
preted as bonding solutions of adatom-CSR couplings. 
This analysis shows also that the adatom-C bonding 
states get successively stronger as the energetical over- 
lap between the TiSR-modified adlevel and the substrate 



CSR's increases. Maximum overlap is achieved for the O 
adatom, for which all bonding adatom-TiSR states lie in 
the lower half of the substrate UVB Q . 

Our new calculations for atomic adsorption on 
TiN(lll) show that the ADOS's [Figs. Ei;b)-(k)] be- 
have in a similar way on this substrate. Figure ^ il- 
lustrates schematically the concerted-coupling model for 
TiN(lll), using the calculated ADOS for fee H on 
TiN(lll), which provides the most complete illustration 
of the different types of adatom-substrate couplings. As 
can be seen, the results for adsorption on TiN(lll) show 
strong similarities to those described above for TiC(lll). 
However, a closer comparison reveals differences: (i) the 
antibonding peak above Ep lies at a lower energy than on 
TiC(lll), +0A eV; (u) the TiSR quenching is stronger 
than on TiC(lll) and involves a larger number of TiSR 
electrons; (iii) the coupling between the TiSR-modified 
adlevel and the XSR's occurs at lower energies than on 
TiC (111). Points (i) and (ii) can be understood from the 
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facts that the TiN(lll) TiSR hes at a lower energy than 
on TiC(lll) and that this causes a larger filling of the 
TiSR on TiN than on TiC. Point (iii) is caused by the 
lower energy of the TiN(lll) UVB, compared to that of 
the TiC(lll) UVB. 

In addition, our state-resolved, real-space, analysis of 
the DOS shows that also on TiN(lll) there is a succes- 
sively increasing dominance of adatom-X bonding states 
[Fig. 121^1)] with decreasing energy. Again, however, dif- 
ferences can be detected: (i) compared to the adatom-Ti 
coupling, the adatom-X coupling is generally stronger on 
TiN(lll) than on TiC (111); (ii) while no adatom-Ti cou- 
pling can be detected in the TiSR energetical region on 
TiC(lll), such adatom-Ti bonding states are found in 
the TiSR energetical region on TiN(lll). Point (i) can 
be understood from the higher ionicity of TiN, which 
results in a stronger LDOS on the N atoms, compared 
to the C LDOS in TiC, and thus a stronger adatom-X 
overlap on TiN(lll) than on TiC(lll). Point (ii) can be 
interpreted as the presence of a weak coupling, similar 
to the adatom-UVB coupling, between the Ti-localized 
CB and the antibonding adatom-TiSR states. Due to 
the CB being partially filled in TiN, this adatom-CB 
coupling results here in filled states, in contrast to TiC, 
where the CB is empty. 

Our results show thus that the atomic chemisorption 
on TiN(lll) can be understood and analyzed in a way 
similar to that on TiC (111). Both types of interac- 
tion of our concerted-coupling model (adatom-TiSR and 
adatom-XSR's) appear to be active on both substrates, 
due to the similarities between the electronic structures 
of the two surfaces. However, some differences are de- 
tected, which can be related to the higher ionicity and 
higher energy of Ep in TiN, compared to TiC. In particu- 
lar, the results indicate that a weak adatom-CB coupling 
is also present. 



IV. DISCUSSION AND CONCLUSIONS 

Our previous analysis shows that, from the DOS per- 
spective, our concerted-coupling model is applicable also 
to adsorption on TiN(lll). In the following, we discuss 
how this model can be used to account for the essential 
features of the calculated trends in ii^ads 0, Hi, which 
varies in a characteristic pyramid-like manner for the 
second-period adatoms B, C, N, O, and F, peaking at 
O, and the third-period adatoms Al, Si, P, S, and CI, 
peaking at S, on both TiC(lll) and TiN(lll) [Fig.E^a)]. 

Given the strong quenching of the TiSR, a strong 
contribution to the bonding can be expected to come 
from the first component of the concerted coupling, the 
adatom-TiSR coupling. Indeed, the adsorption energy 
for O on the TiSR-deficient TiC(OOl) surface is 43% 
weaker than on TiC(lll) [3. The strong bonding contri- 
bution of the TiSR arises from the fact that the antibond- 
ing adatom-TiSR solution lies above Ep, thus remaining 
empty for all adatoms. 



Our calculated ADOS's [Figs. Efb-k)] show that, 
within each adatom period, the TiSR quenching de- 
creases from group-IV to group- VII adatoms, while it 
increases from group III to group IV. These trends could 
indicate an increase in adatom-TiSR coupling strength 
from group III to group IV and a successive decrease from 
group IV to group VII. These trends agree with the cal- 
culated -Eads trends for B ^ C, for Al Si, for O ^ F, 
and for S — > CI. Also, the smaller TiSR quenching in pe- 
riod 3 than in period 2 indicates a weaker adatom-TiSR 
coupling for period 3, in agreement with our -Eads results. 

However, this trend in adatom-TiSR coupling strength 
disagrees strongly with the calculated -Eads trends for C 
^ N ^ O and for Si ^ P — > S. On the other hand, our 
state-resolved DOS's (see Sec. III.B) show that the sec- 
ond component of the concerted coupling, the adatom- 
XSR couplings, gets successively more active as the en- 
ergetical overlap between the adlevel and the substrate 
XSR's increases. As Z increases within each adatom pe- 
riod, the adlevel energy decreases, thus increasing the 
overlap between the adlevel and the energetical region 
of the XSR's (that is, the lower half of the substrate 
UVB). Therefore, the adatom-XSR coupling strengths 
can be expected to increase as Z increases within each 
adatom period, providing an explanation for the overall 
increasing E'ads values from group III to group VI. In ad- 
dition, this provides further explanation for the stronger 
chemisorption of the second-period adatoms, compared 
to the third period, due to the lower free-atom energy, 
within a given group, of the second-period adatoms. 

On the other hand, such a trend is in contrast to the 
-Eads trends for O ^ F and S CI. However, F and 
CI have ionic bonds, as shown by our calculated Bader 
charges, which yield a charge transfer from the surface 
Ti atoms to the F and CI adatoms of 0.75-0.80 electrons, 
and by our calculated electron distributions, which show 
negligible bond charges for F and CI . This implies al- 
most fully occupied outer electron shells for the F and CI 
adatoms, which can therefore be considered to be almost 
inert toward further chemical interaction, that is, toward 
a coupling with the UVB states. Indeed, our ADOS's 
show that the F and CI adlevels do not broaden as sig- 
nificantly as the levels of other adatoms with the same 
amount of overlap with the UVB region [c/., for exam- 
ple, CI and O on TiN(lll)], indicating a weaker adatom- 
UVB coupling. For all the other fee adatoms, the electron 
densities show a clear covalent character for the adatom- 
substrate bond, with the ionic character increasing suc- 
cessively with increasing Z within each adatom period. 
This is confirmed by our calculated Bader charges for 
the adatoms; -1.07, -1.31, -1.34, -1.20, and -0.80 for 
the second-period adatoms B, C, N, O, and F, respec- 
tively (corresponding to 21%, 33%, 45%, 60%, and 80% 
of their empty outer electron shells) and —0.40, —0.89, 
— 1.08, —1.05, and —0.74 for the third-period adatoms 
Al, Si, P, S, and CI, respectively (corresponding to 8%, 
22%, 36%, 52%, and 74% of their empty outer electron 
shells). 
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FIG. 3: Adsorption of the second-period adatoms B, C, N, O, and F and of the third-period adatoms Al, Si, P, S, and CI in 
fee site on TiX(lll). (a) Calculated trends in adsorption energies -Bads (excerpt from Ref. J^). (b)-(k) Calculated densities 
of states (DOS's) for adsorption on TiN(lll). Thin line: adatom-projected DOS. Thick line: difference in total DOS (ADOS) 
before and after adsorption for the surface TiN bilayer, including adatom. Dashed (dot-dashed) line: ADOS projected onto N 
(Ti) atoms of the surface bilayer. Shaded area: energetical location of the substrate UVB. Dashed vertical lines: energetical 
location of the six substrate-UVB peaks. (1) Contour plot of the KS wavefunction at i? = —6.8 eV for the fee O adatom on 
TiN(lll), illustrating a typical bonding adatom-N interaction. 
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FIG. 4: Schematic diagram of the concerted coupling between the adatom state and the TiN(lll) Ti- and N-centered surface 
resonances (TiSR and NSR's), here iUustrated for the case of fee H/TiN(lll). On the left, the DOS for the clean TiN(lll) is 
reproduced. In the middle, our calculated ADOS for H/TiN(lll) is shown, illustrating the quenching and appearance of peaks 
due to adsorption as negative ("— ") and positive ("+") peaks, respectively. Horizontal solid lines represent schematically the 
individual energy levels before and after adsorption. 



Thus, the concerted-coupling model, which well de- 
scribes the calculated DOS results, provides also a way 
for understanding the main features of the calculated 
i^ads trends for both TiX(lll) surfaces. In particular, 
our analysis points toward a greater role for the adatom- 
XSR couplings in this understanding than what has been 
suggested so far in the literature. 

Our results show also a small but interesting difference 
between TiC and TiN: for adatoms with low Z (B and 
Al) chcmisorption is slightly stronger on TiN(lll), while 
for adatoms with higher Z it is stronger on TiC (111). 
For B and Al, the adlevels lie above the UVB re- 
gion. The adatom~XSR contributions to their bonding 
strengths are therefore negligible and their chcmisorp- 
tion is dominated by the adatom-TiSR coupling. Hence, 
our EsAs results imply that the bonding contribution 
from the adatom-TiSR coupling is somewhat stronger 
on TiN(lll) than on TiC(lll), which is in agreement 
with our previous observation that the TiSR quenching 
is stronger on TiN than on TiC. Now, as stated above, 
as the adlevcl energy shifts down and overlaps with the 
UVB region, the adatom-XSR couplings get successively 
stronger. However, for the same adlevel shift, this overlap 
is stronger on TiC than on TiN, due to the smaller UVB- 
CB energy gap of TiC. Therefore, for the same adatom, 
the adatom-XSR contributions to the bonding strength 
are stronger on TiC than on TiN. As Z increases within 



each adatom period, this effect grows strong enough to 

"neutralize" the stronger adatom-TiSR coupling for TiN 
and invert the chemisorption preference between TiN and 
TiC. 

In summary, the application of our concerted-coupling 
model, originally proposed for adsorption on TiC(lll), 
can be generalized to TiN(lll). However, the quantita- 
tive differences in electronic structures between the two 
compounds must be taken into account. This shows that 
adsorption on TiC and TiN is more complex and versa- 
tile than on, for instance, pure metal surfaces and that 
the combination of several different mechanisms must be 
taken into account, giving a potential for variations, very 
useful in applications. 
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